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We present a numerical solution of the four-quark Bethe-Salpeter equa¬ 
tion for a scalar tetraquark. We hnd that the four-body equation dy¬ 
namically generates pseudoscalar poles in the Bethe-Salpeter amplitude. 

The sensitivity to the pion poles leads to a light isoscalar tetraquark mass 
Mcr ^ 400 MeV, which is comparable to that of the ( t // o ( 500 ). The masses 
of its multiplet partners k and ao//o follow a similar pattern, thereby pro¬ 
viding support for the tetraquark interpretation of the light scalar nonet. 
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1. Introduction 

The notion of tetraquarks has undergone a significant paradigm change 
in recent years. There is now increasing evidence for four-quark states in 
the heavy charm and bottom spectrum, with X(3872), y(4260) and the 
charged Z states as the prime candidates [1]. Their internal decomposition 
is under debate: are they ‘compact’ tetraquarks with equally distributed 
constituents? Or do they rather prefer to be in a diquark-antidiquark or a 
meson-molecule configuration? 

The oldest tetraquark candidates are those in the light scalar meson 
nonet: the (t//o(500), «;(800) and ao//o(980), which do not fit into the 
conventional meson spectrum. Why are the oq and /o mass-degenerate, with 
masses close to KK threshold? Why are their decay widths so different from 
a and k? And why are these states so light? In the quark-model counting, 
scalar mesons should be p waves with one unit of orbital angular momentum. 
Hence, they should have masses of about 1.2... 1.5 GeV, similar to their 
axial-vector and tensor counterparts. 

A possible explanation has been proposed long ago [2]: what if these 
‘mesons’ were actually tetraquarks in the form of diquark-antidiquark states? 

* Presented by G. Eichmann at the workshop ’’Excited QCD 2015”, Tatranska Lom- 
nica, Slovakia, March 8-14, 2015. 
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Fig. 1. Four-body BSE for a tetraquark. 


Two scalar diquarks form again a nonet, but the mass ordering is reversed: 
the a would be the lightest state made of up and down quarks only, whereas 
/o and ao would be heaviest and mass-degenerate because they carry two 
strange quarks. This would also explain the decay widths: /o and ao are 
close to KK threshold and therefore narrow, but a and k can simply fall 
apart into two pions or a pion and kaon, respectively, without exchang¬ 
ing gluons. In that case the actual qq ground states would be indeed the 
‘first radially excited nonet’ with masses in the 1.3 ... 1.5 GeV region. The 
non-qq interpretation of the light scalar nonet is supported by a variety of 
approaches such as unitarized ChPT, quark models, or the extended linear 
a model [3]. The status of tetraquarks at large Nc is discussed in Ref. [4]; 
see also Jose Pelaez’ contribution to these proceedings. 

2. Bethe-Salpeter approach 

Irrespective of their internal decomposition, what we can say with cer¬ 
tainty is the following: if four-quark states exist, they should appear in 
QCD’s spectral representation and produce poles in the eight-quark corre¬ 
lator, or more precisely, in the connected qqqq scattering matrix T. Look¬ 
ing for bound-state poles in higher n—point functions is the essence of the 
Bethe-Salpeter approach: at the pole position, the T—matrix becomes the 
product of two Bethe-Salpeter amplitudes, and from the onshell propagator 
pole we can read off the mass of the state. 

The T—matrix can be determined from the scattering equation or inho¬ 
mogeneous Bethe-Salpeter equation (BSE). It has the form T = K+K Gq T, 
where Go is the product of four dressed quark propagators and K is the 
four-quark interaction kernel. The scattering equation at the pole reduces 
to the homogeneous BSE shown in Fig. 1, which has a solution only if 
the T—matrix has a pole. The kernel is the sum of two-, three- and four- 
body irreducible interactions (we already omitted the three- and four-body 
pieces in Fig. 1), and the specific form of the pair interactions is necessary 
to prevent overcounting [5,6]. The equation in the figure is written in the 
‘diquark-antidiquark’ topology (12)(34); there are two further permutations 
(23) (14), (31) (24) with meson-meson topologies. 

The two-body kernel that appears here must be consistent with the 
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underlying quark-gluon structure to preserve QCD’s chiral symmetry. A 
feasible way to achieve this is to employ a rainbow-ladder kernel, where the 
two-quark interaction is simplified to an iterated gluon exchange, and solve 
the Dyson-Schwinger equation (DSE) for the quark propagator with this 
input. As reviewed elsewhere [7], the phenomenology of pseudoscalar and 
vector mesons, as well as baryon octet and decuplet ground states, can be 
reasonably well described in such a setup. This implies not only mass spec¬ 
tra but also their form factors and other properties, and it extends to char- 
monium and bottomonium spectra. Calculations beyond rainbow-ladder 
are underway and aim, for example, to include the ‘pion-cloud’ corrections 
to form factors [8]. 

On the other hand, there are many observables where gluon exchange 
alone does not provide the answer, for instance: scalar and axialvector 
mesons (the ‘p waves’ in the quark model), heavy-light mesons, or excited 
hadrons. When viewed as a qq state, rainbow-ladder produces a it— meson 
mass of 600 ... 700 MeV. Its status beyond rainbow-ladder is presently un¬ 
clear: in some recent calculations the mass goes down even further [9], 
which would again provide support for a qq interpretation. On the other 
hand, these studies do not only produce a light isoscalar but also a mass- 
degenerate light isotriplet, which leads back to the ‘wrong’ mass ordering: 
{a, ao} —)> K —)■ /o. Hence, there are most likely further effects at play. 


3. Tetraquarks: twfo-body equation 

As a first step towards solving the four-body problem, let us take a sim¬ 
pler route that is motivated by the quark-diquark model for baryons [10]. 
When three-body interactions are neglected, the three-body BSE analogous 
to Fig. 1 can be simplified to a quark-diquark BSE where gluons no longer 
appear. Instead, quarks interact with scalar and axialvector diquarks via 
quark exchange. It turns out that the quark-diquark approach, when em¬ 
ploying rainbow-ladder DSE and BSE solutions for its ingredients, provides 
a rather good approximation for the three-body equation regarding nucleon 
and A properties [11]. 

Based on this observation, one can proceed in an analogous way with 
the four-body equation. Fig. 1 can be rewritten as a Faddeev-Yakubovsky 
equation [12], and after assuming separability for the qq and qq scattering 
matrices it reduces to a coupled diquark-antidiquark/meson-meson equa¬ 
tion where the interaction takes place via quark exchange [6]. In fact, the 
equations can be combined into a single meson-meson equation (shown in 
Fig. 2) but not vice versa. This already hints at the meson-molecule nature 
of the system, where diquarks only appear as internal admixtures. Using 
rainbow-ladder DSE and BSE solutions for the ingredients, the solution for 
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Fig. 2. Two-body meson-meson/diquark-antidiquark BSE. Single, double and 
dashed lines are quark, diquark and meson propagators, respectively. The blobs 
are the respective Bethe-Salpeter amplitudes. All ingredients are dressed. 


the a mass was found to be ~ 400 MeV [ 6 ]. 

The evolution of the isoscalar-scalar mass as a function of the input 
current-quark mass is shown in the left panel of Fig. 3. It is noteworthy that 
the omission of the diquark diagram in Fig. 2 changes almost nothing, so the 
diquark admixture is practically irrelevant. The light cj mass is therefore 
tied to the behavior of the light pions which are QCD’s Goldstone bosons in 
the chiral limit: the tetraquark mass lies just below the pseudoscalar-meson 
threshold, except for very low quark masses where it becomes a resonance. 
Moreover, the calculation also predicts a relatively light ssss state as well 
as an all-charm tetraquark below 2m^^. 


4. Tetraquarks: four-body equation 

How reliable are these results given the approximations that were em¬ 
ployed? To find out, let us go back to the original four-body equation. A 
scalar qq state is predominantly in p wave, but a scalar four-quark state 
forms an s—wave orbital ground state. Employing a rainbow-ladder ker¬ 
nel is therefore well motivated from the aforementioned meson and baryon 
studies, and its reliability can be judged from their results. 

The main difficulty in solving the system in Fig. 1 is the complicated 
structure of the tetraquark amplitude. It is a five-point function and de¬ 
pends on four independent momenta. The amplitude can be decomposed 
into two color structures and 256 Dirac-Lorentz tensors, and the result¬ 
ing 512 Lorentz-invariant dressing functions depend on nine independent 
variables. It is practically impossible to solve this equation on present com¬ 
puters, so we must resort to approximations. Since they should respect the 
symmetries of the system, we proceed as follows: 

(i) We keep the 16 (s-wave) Dirac-Lorentz tensors that do not depend on 
any relative momentum. This is a closed set under Fierz transformations, 
so it contains ‘diquark-antidiquark’-like structures such as C ^75 < 8 ) TsC and 
in the (12) (34) topology, but via Fierz identities it also gener¬ 
ates meson-meson-like structures ~ 75 (8 75 in the crossed channels. 
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Fig. 3. Current-quark mass dependence of the scalar-isoscalar tetraquark mass from 
the two-body BSE [6] (left panel) and from the four-body BSE {right panel). The 
dashed lines are the respective thresholds. 


(ii) We take into account both color structures, which can be written in 
either of the Fierz-equivalent topologies: diquark-antidiquark (3 ( 8 ) 3, 6 ( 8 ) 6 ) 
or meson-meson (1 ( 8 ) 1 , 8 ( 8 > 8 ). 

(iii) We arrange the nine momentum variables into irreducible multiplets 
of the permutation group ^4 [13], which allows us to switch on groups of 
variables separately in the solution process and thereby judge their impor¬ 
tance. The nine variables can be grouped into a symmetric singlet 5o, a 
doublet Vq, and two triplets 71, Ti: 


‘^0 — ^ {p^ + + k^) , V — 


45n 


\/3 (g^ - P^) 

— 2k‘^ ’ 


( 1 ) 


where g^, k"^ are the Mandelstam variables for the three topologies. 
Hence, the dressing functions can be written as /* = /i(5o, 71, Ti)- The 

doublet phase space that remains invariant under the equation forms the 
interior of a triangle bound by the lines = 0 , g^ = 0 and = 0 (see left 
panel in Fig. 4), whereas the triplets form a tetrahedron and a sphere. 

The resulting eigenvalue of the four-body BSE is plotted in Fig. 4. Since 
the equation has the form = AT, one can extract the mass of the state 
by tracking the largest eigenvalue of the kernel as a function of the total 
momentum squared: A(P^ = —M^) = 1. If we allow the system to generate 
a momentum dependence in 5o only, the eigenvalue (second curve from 
below) would cross the line A = 1 at M ~ 1400 MeV. This is essentially 
‘four times the constituent-quark mass’ and would describe a rather simple 
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Fig. 4. Left panel: Mandelstam triangle whose axes (not shown) are the doublet 
variables in Eq. (1). The contour plot refers to the magnitude of the leading dressing 
function. Right panel: BSE eigenvalue after switching on different multiplets. 


tetraquark. Switching on the triplets does not change the result appreciably, 
but adding the doublet V has a drastic effect: the eigenvalue curve is now 
almost flat and crosses at M ~ 400 ... 500 MeV. 

What happened? The exterior of the triangle exhibits pion and diquark 
poles at timelike values of the Mandelstam variables, illustrated by the 
dashed and double lines in Fig. 4. We emphasize that the four-body equation 
knows a priori nothing about pions and diquarks but rather generates them 
dynamically. They influence the behavior of the dressing functions in the 
interior, which is clearly visible in the contour plot: the pion mass is small 
and the pion poles are close to the triangle, whereas the scalar diquarks 
have no visible effect because their mass scale ~ 800 MeV is much larger. 
If the pion poles are removed by hand via phase space cuts, the uppermost 
eigenvalue curve in the right panel approaches again its companions. 

Therefore, we can confirm our findings from the two-body equation: the 
system is dominated by pions rather than diquarks. While this would be 
difficult to judge from the tensor structures alone (because they mix under 
Fierz transformations), the momentum dependence of the amplitudes clearly 
provides such information. The equation also produces a physical threshold: 
if M > 2 m 7 r, the poles enter the integration domain and the tetraquark 
becomes a resonance. In that case we must resort to extrapolations to 
estimate the real part of the mass. The current-mass dependence of the 
scalar-isoscalar tetraquark mass is shown in the right panel of Fig. 3 and 
mirrors the result from the two-body equation. At the physical u/d mass, 
our results for cr, k and ao//o are 


Mcr ~ 380 MeV, 


Mk ~ 700 MeV, 


M„,/;„~ 920 MeV. (2) 
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5. Summary 

We have solved the four-quark Bethe-Salpeter equation for a scalar 
tetraquark. The outcome provides a simple dynamical explanation for the 
light a mass: the system is dominated by pion poles, which are dynami¬ 
cally generated in the solution process and drive the tetraquark mass from 
‘four times the constituent-quark mass’ to M„ ~ 400 MeV. The masses for 
the multiplet partners k and ao//o follow an analogous pattern. This pro¬ 
vides support for the interpretation of the light scalar nonet as tetraquarks. 
However, they are not diquark-antidiquark states but predominantly ‘meson 
molecules’, which explains their mass ordering as an indirect consequence 
of QCD’s spontaneous chiral symmetry breaking. 
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